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We present an efficient method to solve the quantum discord of two-qubit X states exactly. A 
geometric picture is used to clarify whether and when the general POVM measurement is superior 
to von Neumann measurement. We show that either the von Neumann measurement or the three- 
element POVM measurement is optimal, and more interestingly, in the latter case the components 
of the postmeasurement ensemble are invariant for a class of states. 
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The notation of quantum discord was proposed in 2001 [1] . It is regarded as a measure of quantumness of correlation, 
even in the absence of quantum entanglement. For ten years, many works have been devoted to the significance and 
application of quantum discord (see, for example, [2]). The analytical expressions for quantum discord have been 
obtained only in a few cases including two-qubit Bell-diagonal states [3], rank-two states [4] and Gaussian state [•")]. 
However there is no exact result so far for two-qubit X states (i.e., the states such that the non-zero elements of the 
density matrix only lie along the diagonal or skew diagonal). In this paper we present an efficient method to solve 
this problem. 

As we have known, the major difficulty in the calculation quantum discord is how to acquire the maximal information 
about one particle by measuring the other particle. Given a bipartite state p^^ , perform on particle A a generalized 
measurement, denoted by POVM M = {Mi} with Af; ^ J^i ^'^i — 1- The accessible information about particle B 
is given by S{p^) — S{p^\M). Here the conditional entropy S{p^\M) = ^^^PiS^p^^^') is the weighted average of 
the states = Tr^[(Mi ig) l)p^^]/pi that correspond to the individual outcomes with probabilities pi = Tr[{Mi (g) 

l)p'^^] — Tr{A'Iip^). The maximization over all POVMs gives the classical correlation, C = maxM [Sip^) — S{p^\M)'\ . 
Quantum discord Q is given hy Q ~ X ~ C, where I is the total correlation quantified by the mutual information, 
I — S{p^) + S{p^) — S{p^^). It is a formidable task to find the optimal measurement among all M to achieve 
the minimal value of the conditional entropy S{p'^\M). Much effort [()-9], analytical or numerical, has been made 
in studying the optimization for two-qubit states. However there is no definite answer as to whether and how the 
quantum discord is determined by the general POVM measurements. 

The measurement M on A induces the decomposition of p^ into the ensemble {pi, p^'*^'}. For two-qubit states, we 
have known that all are distributed, in terms of the Bloch vectors, in an ellipsoidal region in three-dimensional 

real space. This region is called quantum steering ellipsoid [1(1], which we denote by £. It has been shown that this 
geometric picture is very useful in the discussion of the quantum discord of two-qubit states [4, 11]: We need only 
consider the decomposition with the form of p^ = J2iPiPf^ where all pf are distributed on the surface of £. The 

problem is then to find the minimal value of the average entropy = ^iPiS{pf), which we denoted by <S'^i„. As we 
show later, we benefit greatly from this geometric picture in the case of two-qubit X states: The optimal measurement 
on A, or the optimal decomposition of p^ , can be determined unambiguously, and thus the exact result of quantum 
discord is attained. 

Note that there are infinitely many states corresponding to a given £. Denote by [p"^'^]^ the set of all X states 
having the identical <B. We show that all steering ellipsoids associated with X states are classified into three types, (i) 
€.(_>.: For any state in the set [jO'*^]e^, the optimal measurement is such a von Neumann measurement that induces a 
horizontal decomposition of p^ in the geometric picture, (ii) €j: The optimal (von Neumann) measurement induces a 
vertical decomposition of p^ for any p^^ S [p'^'^](ej. (iii) £a: For some states in [p'^^Je^, the optimal decomposition 
is horizontal, while for others the three-state decomposition is optimal. More interestingly, for the different states in 
the latter case, say p"^^ and p'^^ , the optimal decomposition is realized on the identical components, that is, the 
reduced states p^ and p'^ are decomposed optimally as Y^^=iPiPF ^^'^ Y^'i=iP'iPF^ respectively. 

In other words, the optimal decomposition is mainly determined by the property of £. Although there is no a 
systematic way to describe the optimal measurement, the optimal decomposition does be described clearly in the 
geometric picture. 
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Any X state, up to local unitary operations, can be written as 



Vu 

where u,v ^ and satisfy ^ ad and v"^ ^ be. The steering eUipsoid € is given by + + ~ ^0)^/^3 = 1, 

where £1 = "+•" ^ £„ — |m-^;| ^ ^ _ , ^"^w^^L and zn = ^-ffw^r^- The reduced state corresponds to 

^(a+&)(c+d) ' ^(a+f,)(c+d) ' (a+b)(c+d) U (a+6)(c+rf) 

the point B = (0, 0, a — 6 + c — d), which is on the z axis and in the interior of £. Note that if u = u or ad = 6c the 
ellipsoid degenerates to an ellipse or a line segment. We claim in advance that the procedure presented below can be 
readily applied to these degenerate cases. 

For the state given by (1), the optimal von Neumann measurement must lie in the x-z plane [.'f]. It is not difficult 
to see that the same conclusion holds for POVM measurements: Each operation element of the optimal POVM 
measurement must lie in the x-z plane. It means that we need only analyze the decomposition of into the convex 
combination of the stats on the ellipse given by x^/^f + (z — zo)^/^3 — Ij which we denote by E. All available 
decompositions can be described geometrically as follows. In the ellipse E, we plot an inscribed polygon B1B2 ■ ■ ■ Bn 
encompassing the point B (Fig. 1(a)). Then decomposing p^ into ^^Pipf amounts to expressing point B as the 
convex combination of the points Bi, B2, ■ ■ ■ , Bn- This case corresponds to the n-element rank-1 POVM measurement 
performed on A. In fact, we need only consider the case of n ^ 4 [12]. Any chord of the ellipse passing through point 
B denotes a two-state decomposition which comes from von Neumann measurement on A (see Fig. 1(a)). 



/a u\ 
b V 




FIG. 1. Schematic plot of ellipse E. Point B denotes the reduced state . In (a), the line segment T1T2 denotes an tilted 
two-state decomposition of . The inscribed quadrilateral B1B2B3B4 denotes a four-state decomposition of p^ . In (b), the 
horizontal line segment EF and the vertical GH indicate the horizontal and the vertical decomposition of p^ , respectively. 

In the case that only von Neumann measurements are allowed, Ali et al [(i] claimed that for the X state (1) 
the optimal observable is either or cr^, which give rise to the horizontal or vertical decomposition, respectively 
(Fig. 1(b)). Denote by Sx the von Neumann entropy of the state corresponding to the point X. The horizontal 

decomposition gives the average entropy 5^ = peSe + PfSf = Se = Sp, and the vertical one gives = pgSg + 
PhSh- All's conclusion is — minlS"^, S^}. However, it is pointed out that there exist some states for which 
neither nor is optimal [7]. A tilted decomposition, which comes from measuring the observable such as 
(j^cosO + afsinO^ will give a smaller average entropy. But it is not the end of the story. Is the two-state tilted 
decomposition optimal? What about three- or four-state decomposition? These problems remain open. 

To solve these problems, let us point out an evident fact. The tilted decomposition of is equivalent to a trivial 
four-state decomposition. In fact, we can plot another tilted line r3T4 (obtained by rotating 180° around the z 
axis) which gives rise to the same average entropy as T1T2 (see Fig. 2(a)). Therefore the four-state decomposition 
{Ti, T2, T3, T4} has the same effect as {Ti,r2}- This observation leads us to the general four-state decomposition 
as illustrated in dashed lines in Fig. 2(a). Let B1B3B2B4 be an isosceles trapezoid inscribed in the ellipse with 
BiBi and B2B3 parallel to x axis. Point B is in the interior of the trapezoid. Then the average entropy is given by 
S'4 = X]i=i Pi'S'si • Point B may not be the intersection point of the diagonal lines of the trapezoid. 
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FIG. 2. In (a), the four-state decomposition {Ti}i=i,... ,4 is equivalent to the 2-state one {Ti,T2}, while {-Bi}i=i, -. ,4 is more 
general. In (b), we plot schematically the curve S^{z). Two tangent points indicate the solution of (2). 



Define the function h{x) as h{x) ~ log2 — log2 for x e [0, 1]. We have Sb^ — h{OBi), where OB^ 

is the distance from the origin point O to point Bi. Since OBi = OB^, OB2 — OB3, pi ~ and p2 ~ P3, we rewrite 
^4 as 5*4 = 2piSbi + 2p2*S'_B2- Let zb, z\ and Z2 denote, respectively, the z-coordinate of B, B\ and i?2. Note that 
Sbx can be expressed as S'<^(zi), which is the average entropy given by the horizontal decomposition of the state 
corresponding to point (0, zi) in Fig. 2(a). Similarly, = 6*^(22). With the probabilities given by p^ = of^^"^ ^ 

and P2 ~ N , the problem is to minimize ^4 = 2pi5^(zi) + 2p2'S'^(z2) over zi £ \zh^ zb\ and Z2 € [zb, -zg] with 

zh and zq the coordinates of the lower and the upper vertex of the ellipse E respectively. Taking partial derivatives 
of 5*4 with respect to zi and Z2, we have 

dS^jzi) _ dS^{z2) _ S^{zi) - S^{z2) , , 

dzi dz2 Z2 

It means that if there exists a four-state optimal decomposition, then the function S^(z) has at least two inflection 
points, that is, the equation S ^{z) / dz"^ = has at least two solutions (see Fig. 2(b)). This observation motivates us 

to investigate the properties of the horizontal average entropy S^{z) — h{r{z)) with r(z) = [z^+^4[l — (z — zo)^/^§] 
for z € [z/f , zg\- We have the following lemma. 

Lemma. The horizontal average entropy Sm^{z) has at most one inflection point. 

This lemma can be proved by directly analyzing the properties of the second derivative <P S ^^{z) / dz"^ . But the 
proof is too technical to be described here. Instead let us consider a concrete example: a parameterized X state. Let 
p"^^{ki, ^2) is such an X state that comes from (1) with u and v replaced by ki\/ad and k2\fac respectively. Assume 
that fci, ^2 G [0, 1]. When ki — k2 = 1, the ellipse E is inscribed to the unit circle. It is the largest ellipse when a, b, 
c and d are fixed. With fci and /c2 decreasing from 1 to 0, the ellipse shrinks to the z axis while the upper and the 
lower vertex remain unchanged. Note that if ki — k2, this process is similar to that of an X state undergoing local 
dephasing channel. Now we let a = 0.6717, b — c = 0.125 and d = 0.0783. When fci = and fc2 = 0.8, the state, up to 
local flip operations, is just the one that has been used in [7] to show that the tilted decomposition is superior to both 
horizontal and vertical decomposition. In the following we let k\ = k2 = k for simplicity. For the sake of explicitness, 
we consider the difference A(z, k) = S^{z, k) — S^{z, k). As S^{z, k) is a Hnear function of z, A(z, k) has the same 
convexity as S^{z, k). In Fig. 3, we plot A(z, k) for five different values of k. With k decreasing, we see that A(2, k) 
transforms continuously from the convex to the concave, and meanwhile there is at most one inflection point. The 
same is true for Si^{z). 

Now according to Lemma, Eq. (2) cannot be satisfied. 

Proposition 1. For all X states, the four-state decomposition can not be optimal. 

The remaining cases are two- and three-state decomposition. Both of them are regarded as the extremal cases of 
four-state decomposition. Two-state decomposition corresponds to (i) zi = Z2 = zb (horizontal), or (ii) zi = zu and 
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FIG. 3. (Color online) Plot of A(2,fc) with k taking the values (from bottom to top) of 0.2839, 0.2827, 0.2822, 0.2817 and 
0.2805. The bottom curve (red) is convex, and the top (blue) is concave. The three middle curve (black) are neither convex 
nor concave. Each of them has a unique inflection point (denoted by the dot). 



Z2 — zq (vertical). Three-state decomposition corresponds to (iii) Z2 — zq and zi G {zh,zb) (upper-triangle), or (iv) 
Zi = zh and Z2 € {zbt^g) (lower-triangle). If S'^(z) is convex or concave, it is easy to see that case (iii) and (iv) 
cannot satisfy (2). Then we have 

Proposition 2. // the horizontal average entropy S^{z) is a convex (concave) function for z G \z}i,zq\, then the 
ellipsoid is the type of (<B^): for all X states associated with the ellipsoid, the optimal decomposition is the 
horizontal (vertical) one. 

In proving the Lemma, we see that (i) if u + v ^ \ad—bc\, S^{z) is a convex function, (ii) if {u + v)"^ < (a — 6)(d— c), 
S^{z) is a concave function, and (iii) if the center of the ellipse coincides with the origin (i.e., zq — 0), then S^{z) 
is either convex or concave. Then from Proposition 2 we have two corollaries. 

Corollery 1. If u + v ^ \ad — bc\, the horizontal decomposition is optimal. If {u + w)^ ^ {a — b){d — c) , the vertical 
decomposition is optimal. 

Corollery 2. // the center of the ellipse coincides with the origin, then the horizontal decomposition is optimal for 
oblate ellipse (i.e., £i > i^), and the vertical decomposition is optimal for prolate ellipse (i.e., £i < £3). 

Note that the same conclusion as Corollary 1 has been obtained in [!)], and that Corollary 2 is not restricted to the 
states with maximally mixed marginals. 

Something interesting is going on here. Let us consider the case that S^{z) has a unique inflection point. The 
pictures of such a sort of S^{z) are plotted schematically in Fig. 4(a)-(c), which in fact correspond to the three 
middle curves in Fig. 3. In the cases illustrated in Fig. 4(a) and (b), we can always draw a line, denoted by Lq, 
tangent to the curve S^{z) and passing through its right end point. Denote by z* the z-coordinate of the tangent 
point. We see that z* satisfies 



That is, z* is the solution of Eq. (2) in the extremal case of Z2 — zq. If zb G {z* , zq), this solution gives rise to the 

three-state decomposition of , which leads to the minimal average entropy ^jj^jj^. Fig. 4(d) illustrates the upper- 
triangle decomposition: the horizontal line PQ denotes the horizontal decomposition of the point K — (0, z*), and the 

optimal decomposition of B is given by {G, P, Q}. It follows that = p* S^{z*) + {1—p*)Sg — P* Sp + {1—p*)Sg 
with p* = f^^f?-. We emphasize the fact that the solution z* is independent oi zb- It implies that for any point 
B G KG, the optimal decomposition is always given by {G,P,Q}. Of course, for B' S HK, we cannot benefit from 
the three-state decomposition and the optimal choice for B' is still the horizontal one. 

In the cases illustrated in Fig. 4(b) and (c), we can also draw a tangent line Ljf from the left end point. The 
z-coordinate of the tangent point is the solution of (2) in the extremal case of zi = zh. This solution also provides a 



dSM 



S^jz*) - S^jza) 

Z* - ZG 



(3) 



dz 
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FIG. 4. Fig. (a)-(c) are the schematic plots of S^{z) with only one inflection point (solid curve). S'j(z) is labeled by dotted 
line. Fig. (d) illustrates the optimal upper-triangle decomposition {G,P,Q}. 



three-state decomposition. But it cannot be optimal, because the average entropy S given by this decomposition is 
larger than that given by the horizontal and vertical one. 

It should be mentioned that the above analysis, in particular Fig. 3 and Fig. 4, is sensible in the case that the 
center of the ellipse is on the upper half of z axis, namely, zq > 0. If zq < 0, analogous analysis shows that in this 
case we will encounter the lower-triangle decomposition. From these considerations, we have 

Proposition 3. If S^{z) is not larger than S^{z) for all z e {zh,zg), and if S^^{z) is neither convex nor con- 
cave, then the ellipsoid is the type of 2;a - some X states associated with the ellipsoid admit the optimal three-state 
decomposition. Moreover, the components of the optimal three-state decomposition are the same for these X states. 

Let us discuss some applications of our results. Consider first the X states with the property that sf^ — . 
This sort of state appears on the boundary of the maximally discordant mixed states (MDMS) [l-H], and also on the 
boundaries of the relationship between the discord and the mixedness measured by the von Neumann entropy or linear 
entropy [14]. Proposition 3 shows that for these states the optimal measurement must be three-element POVM. 

Another scenario in which our results are relevant is that of the sudden transition between classical and quantum 
decoherence presented in [l-l]. It has been shown that the two-sided phase damping channel will cause the steering 
ellipsoid shrinking towards the z axis [11]. The ellipse E transforms from an oblate ellipse (if exists initially) to a 
circle, then to a prolate ellipse. Note that the states used in [l-j] to demonstrate the sudden transition are the X states 
with the maximally mixed marginals, which fall under the class of states with the property of zo = 0. According to 
Corollary 2, the sudden transition follows from the sudden change of the optimal decomposition from the horizontal 
to the vertical. However, for general X states with zq ^ 0, the optimal decomposition changes continuously, rather 
than suddenly, from the horizontal to the vertical. We think that the sudden transition cannot appear in this case. 

In summary, we present an efficient method to solve the quantum discord of two-qubit X states, and give a conclusive 
answer as to what type of and under what conditions the POVM measurement is optimal. More interestingly, for a 
class of states the postmeasurement ensembles of qubit B induced by the optimal three-element POVM measurement 
on A have the identical components. The geometric picture developed in this work is useful in discussing the boundary 
states such as MDMS and the dynamics of quantum discord. The remaining question is, under what conditions is the 
ellipsoid the type of or £a? Also, can this method be generalized to more general two-qubit states? These 

questions will be discussed elsewhere. 



6 



This work is supported by National Nature Science Foundation of China 10875053, the CAS, and the National 
Fundamental Research Program 2007CB925200. 



[1] H. OUivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 (2001); L. Henderson and V. Vedral, J. Phys. A: Math. Gen. 
34, 6899 (2001). 

[2] V. Vedral, Phys. Rev. Lett. 90, 050401 (2003); A. Datta, A. Shaji, and C. M. Caves, Phys. Rev. Lett. 100, 050502 
(2008); A. Shabani and D. A. Lidar, Phys. Rev. Lett. 102, 100402 (2009); K. Modi, T. Paterek, W. Son, V. Vedral, 
and M. WiUiamson, Phys. Rev. Lett. 104, 080501 (2010); D. Cavalcanti, L. Aolita, S. Boixo, K. Modi, M. Piani, and 
A. Winter, Pirys. Rev. A 83, 032324 (2011); V. Madhok and A. Datta, Phys. Rev. A 83, 032323 (2011); A. Streltsov, 
H. Kampermann, and D. Brufi, Phys. Rev. Lett. 106, 160401 (2011); M. Piani, S. Gharibian, G. Adesso, J. Calsamiglia, 
P. Horodecki, and A. Winter, Phys. Rev. Lett. 106, 220403 (2011); M. F. Cornelio, M. C. de Oliveira, and F. F. Fanchini, 
Phys. Rev. Lett. 107, 020502 (2011). 

[3] S. Luo, Phys. Rev. A 77, 042303 (2008). 

[4] L.-X. Cen, X.-Q. Li, J. Shao, and Y. Yan, Phys. Rev. A 83, 054101 (2011); M. Shi, W. Yang, F. Jiang, and J. Du, ,J. 
Phys. A 44, 415304 (2011). 

[5] P. Giorda and M. G. A. Paris, Phys. Rev. Lett. 105, 020503 (2010); G. Adesso and A. Datta, Phys. Rev. Lett. 105, 030501 
(2010). 

[6] M. Ali, A. R. P. Ran, and G. Alber, Phys. Rev. A 81, 042105 (2010). 
[7] X.-M. Lu, J. Ma, Z. Xi, and X. Wang, Phys. Rev. A 83, 012327 (2011). 
[8] D. Girolami and G. Adesso, Phys. Rev. A 83, 052108 (2011). 

[9] Q. Chen, C. Zhang, S. Yu, X. X. Yi, and C. H. Oh, Phys. Rev. A 84, 042313 (2011). 
[10] F. Verstraete, Ph.D. thesis, Katholieke Universiteit (2002). 

[11] M. Shi, F. Jiang, C. Sun, and J. Du, New Journal of Physics 13, 073016 (2011). 
[12] G. M. D'Ariano, P. L. Presti, and P. Perinotti, J. Phys. A 38, 5979 (2005). 
[13] F. Calve, G. L. Giorgi, and R. Zambrini, Phys. Rev. A 83, 012102 (2011). 

[14] A. ALQasimi and D. F. V. James, Phys. Rev. A 83, 032101 (2011); D. Girolami, M. Paternostro, and G. Adesso, J. Phys. 
A 44, 352002 (2011). 

[15] J. Maziero, L. C. Celeri, R. M. Serra, and V. Vedral, Phys. Rev. A 80, 044102 (2009); L. Mazzola, J. Piilo, and 
S. Maniscalco, Phys. Rev. Lett. 104, 200401 (2010). 



